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Abstract 

We consider a thought experiment, in which a neutrino is produced by an elec- 
tron on a nucleus in a crystal. The wave function of the oscillating neutrino is 
calculated assuming that the electron is described by a wave packet. If the electron 
is relativistic and the spatial size of its wave packet is much larger than the size 
of the crystal cell, then the wave packet of the produced neutrino has essentially 
the same size as the wave packet of the electron. We investigate the suppression 
of neutrino oscillations at large distances caused by two mechanisms: 1) spatial 
separation of wave packets corresponding to different neutrino masses; 2) neutrino 
energy dispersion for given neutrino mass eigenstates. We resolve contributions of 
these two mechanisms. 



1 Introduction. 

There are two different approaches to neutrino oscillations in the literature: one of them 
deals with the wave function of free neutrinos, while the other considers the propagator 
of virtual neutrinos. The latter approach was analyzed in the papers [I], [2j for a thought 
experiment, where a neutrino was considered to be produced by an electron on the target 
nucleus A, and captured by the nucleus B in the detector. The process was described as a 
two-stage Feynman amplitude, e + A — > C + u, v + B — > D + 1 with a virtual neutrino, its 
Green function connecting the production and detection points: xa = (*a, x .a) and xb = 
(£b, xb). As it is well known, at a large distance Tab = |xb— from the production point 
virtual particles become effectively real and one may speak about their wave function, in 
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particular, about the neutrino wave function, ip v {x). Though this statement is well known 
and practically evident, an explicit expression for ip v can be instructive for the description 
of the effect of oscillation suppression in a thought experiment. 

There are two mechanisms of erasing oscillations. The first one is spatial separation 
of neutrino wave packets of different mass eigenstates (see refs.^j [Hj). The second 
mechanism is caused by neutrino energy dispersion (see e.g. ref.jZj). The amplitude of 
the two-stage process e + A — > C + Vi and z/j + B —> D + 1 with a virtual neutrino of given 
mass rrii is determined by the standard rules of quantum field theory pQ : 

At = J d 4 x 1 d 4 X2ipt{x 2 )^l ) {x2)ipB{x2)Gi(x2 - x 1 )^c{x 1 )^ e {x 1 )^ A {xi), (1) 

where x\p are 4-dimensional coordinates, Xk = (£&,Xfc), Gi(x 2 — x\) is the Green function 
of the i-th neutrino mass eigenstate, and non-essential spin factors are neglected. The 
term = Ai A* in the probability was integrated over the phase space of the final 
particles in ref.pQ. After integration the interference term with i ^ j vanishes at large 
spatial separation |x^ — x#|. 

In this approach, however, we were unable to resolve the contributions of the two 
mechanisms. In this note we try to separate the "Siamese twins". For this purpose we 
consider only the first stage of the process, e + A — > C + u, with a free neutrino. We 
consider the case when the spatial size of the electron wave packet is much larger than 
the size of the crystal cell, which determines the localization of the nucleus. The opposite 
case will be described elsewhere. 

2 Wave function vs amplitude 

The non-normalized wave function of a neutrino produced in the reaction e + A — > C + v 
is 

i 

ipi(t, x) = J d A xiGi{x - xi)4)* c {xi)4) e {xi)4) A {xi), (2) 

where x = (t, x) is the space-time coordinate and \Pj) is the i— th neutrino mass eigenstate. 

In this equation the product 'ipc( x i)' l Pe{ x i)' l l ; A(xi) serves as a local source of neutrinos. 
To calculate the amplitude of neutrino interaction with the nucleus B one would evidently 
substitute this expression for ipi into the integral over x with the wave functions of other 
particles participating in the reaction, according to eq. (|T]) . This naturally agrees with the 
general prescription. 

Like in our recent paper pP we consider the initial nucleus bound in a crystal and 
describe it by a stationary wave function localized near the point x = : 

* A (x) = F A (x) e~ itEA , (3) 
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where Ea is the energy of the nucleus. The Fourier transform of i^(x), which is required 
in what follows, is 

K A {<U) = J rfxF A (x)e- iq - 4X . (4) 

By assumption, the nucleus A is at rest and, thus, Ka(<Ia) is centered at = with 
the uncertainty a a ~ a -1 - 

The wave function of the incident electron is considered to be a wave packet: 

* 6 (ar) = J dq e K e (q e - Pe ) e *fe(*-*«)-iE«(qe)« = e ; Pe (x-x e )-iE e ( Pe )i ^ ( x _ e ^ _ ^ 

Here e = p e /p e , p e = |p e |, E e (q e ) = \/ql + mf, the Fourier amplitude i^ e (q e — p e ) 
is centered near q e = p e with the uncertainty a e , the center of the packet envelope 
F e (x — ev e t) is at the point x e at the moment t = 0, and the electron group velocity t> e 
is defined as 

= <9E e (q e ) _ Pe r ^ l _< f gN 

In what follows we assume the electron to hit the nucleus A at t = and the collision to 
be central, i.e. x e = 0. 

The recoil nucleus C is described by the plane wave 

unless its momentum is comparable with a a, which is an extremely rare case. 
For the Green function the following expression can be derived 

1 /*oo 

G;(t,x) = — / ^ e -^+V- 2 - m ?M. (8) 

4tt|x| 

Let us now substitute expressions ©, ©, (J7J) an d (JHJ) into eq.(J2J) and perform trivial 
integration over x\\ 

ipi(t, x) = i J dudq e K e (q e - p e ) dq A K A (q A ) 

5(E e (q e ) + E A - E c - uj)5(q e + q A - p c - k;) exp(^r - iut) . (9) 
Here r = |x|, kj = nfej, n = x/r, hi = sJoS 2 — mf and the expansion 

|x — xi| ~ r — nxi, (10) 

is performed assuming r to be much larger than the interaction region. In eq.tjHJ) and in 
what follows we omit some non-essential numerical factors. We integrate eq.© over q^ 
and L) and obtain 



i>i(x) 



^ J dq e K e (q e - p e ) K A (— q e + Pc + kj) exp [iA^r - iwt] , (11) 



3 



where now u;(q e ) = E e (q e ) + E A — E c . 

While proceeding with the calculations we bear in mind the range and hierarchy of 
the quantities involved: 

rrii < cr e < o A < m e < p e < M A , M c . (12) 
In particular, we consider the case when Ka{q.a) is wide in comparison with K e (q e ): 

a e < o A - (13) 

This allows to set K A (— q e + Pc + kj) = (— Pe + Pc + k°), k° = kj | qe=Pe throughout 
the essential range of integration over q e . Assuming that the momentum distribution of 
the electron is sufficiently narrow we may expand the integrand in terms of q e near the 
central electron momentum p e : 

u;(q e ) = u° + f e e(q e - p e ), 

k i (q e ) = k? + —e{q e -p e ), (14) 
where we introduce the neutrino group velocities analogously to the electron one: 



the upper index "0" means that the corresponding quantities are calculated at q e = p e . 

Taking into account eqs.tjH}, (fTTj) and ()14j) we obtain the following simple expression 
for the wave packet of the produced neutrino through the envelope of the wave packet of 
the incoming electron in coordinate space, F e : 

e ik°r-iiu t / y \ 

^(t,r,n) = K A {-p e + p c + k?n)F e I ^■(r-v i t)ej . (16) 

The factor v e /v,i ~ 1 makes the neutrino wave packet a little bit wider then the electron 
one. It is not essential for our purposes and will be omitted in what follows. 

Equation . which is one of the main results of this paper, is quite natural. If there 
were a long wave packet of the incoming electron, it would create a packet of neutrinos 
with a similar length. A good analogy is the scattering of a sound wave on a target which 
creates another sound wave. The duration and, correspondingly, the spatial length of the 
produced wave packet should be equal to the duration and size of the original one. 

Strictly speaking, one has to calculate the amplitude ((H) to determine the probability 
of the oscillating behavior of neutrinos. However, one may rely on a simplified approach 
based on the interpretation of the absolute value squared of the neutrino wave function 
as the probability density for the particle to be found at a spatial point x at a given 
time t. Such an approach is valid, if we deal with wave packets, the longitudinal size of 
which is much larger then their wave length, and their transversal size is much larger than 
their Compton wave length (see, e.g., [Hj, [01)- If we neither register the nucleus C nor 
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measure the time of neutrino detection, we are interested in the detection probability of 
the neutrino v\ at point x: 



P 



(x) 



dpcdt 



(17) 



For simplicity, in what follows we assume i — 1,2, / = e, /i, U e \ = cos 8, U e2 = sin 9. 
For the production probability we obtain 

Pu^vM = ^(sm20) 2 J dp c K 2 A {-p e + pc* + fc?n) 
y dtF e ((r - ^)e) F e ((r - « 2 t)e 



l- 7 cos 



r<5m 2 



2cj° 



(18) 



wij — mj. 



where / = J dtF 2 ((r — t)e), 5m 2 

The term proportional to cos(r<5m 2 /2ci; ) in eq.(|18|) describes oscillations. It vanishes 
at large distances for two different reasons. 
1. Packet separation. 

When (t> 2 — V\)r > l/cr e , that is 



I J-'osc ) J-'osc r- r> 

o P om z 



(19) 



the product F e (r — vit)F e (r — V2t) is nearly zero for every t, and the oscillating term 
vanishes even before integration over pc- 
2. Neutrino energy dispersion. 
Note that u° depends on p c : 



E e {p e ) + E A -M C - 



P 2 c 



2M, 



c 



(20) 



where M denotes the C nucleus mass. Owing to the factor K\{— p e +pc+k^n) in eq. (|T%|) . 
the effective size of the region of integration over pc is of the order of a a- Thus u° varies 
near its central value (which is roughly equal to p e ), and the variation equals p^/Mc. 
Here, the central value p^ is determined by the equation 



Fe 2p e 2M C 



n. 



(21) 



If the angle between vectors e and n is sufficiently large (|e — n\p e 3> Mq — Ea, p 2 /2Mc) 
then 



p c ~ p e \e - n| 



(22) 

This is the case considered below. The variation of uj° equals p e |e — n\a A /Mc- Thus 
integration of cos (rSm 2 /2u>°) over p<^ in eq. ()18[) leads to the vanishing result if 



We see that there are two competitive mechanisms for suppression of neutrino oscil- 
lations. If 

» e |e — n| , A . 

v e < cr A m ' , (24) 
Me 

then the energy dispersion mechanism dominates, and 

j^ = L *>. (25) 
e — n (7a 



If 



aA ~T7 < °"e < cr A , (26) 



then the packet separation works, and 

T - T — (97) 

-'-'sup - LJ osc ■ V / 

These results coincide with those obtained in ref. pQ. 

One more comment is worth making at this stage. As we said above, the interference 
disappears when two neutrino wave packets ipi an d if)j, eq.f|16|) . cease to overlap. At first 
sight this statement is at odds with the expression for the production amplitude (JTJ) of 
lepton I on the nucleus B. Indeed, the product of amplitudes A4 and A*j, which enters 
the probability of the process (see eqs.(27)-(29) of ref. P), does not vanish even when the 
product %l)i(x)%l)*(x) vanishes, because the amplitude contains an integral over x and the 
product of integrals never vanishes. However, one can check that after integration over 
the phase space of the final particles the product of the integrals vanishes exactly when 
the neutrino wave packets no longer overlap. Such integration over the phase space makes 
the result effectively local. 

This conclusion is intuitively clear for the following reasons. The product of the 
amplitudes prior to integration over the phase space describes the production probability 
of a plane wave final state, because the final states are taken as momentum eigenfunctions. 
It is evident that such a probability never vanishes even in the case of infinite separation 
of neutrino wave packets. It is essentially the same as the excitation of a resonator by two 
wave packets. A resonator with a very large Q-factor would stop to oscillate only after 
a very long time. So, if such a resonator is hit by one wave packet and after a while by 
another delayed wave packet, the interference between the packets would still be observed 
by such a resonator because it keeps oscillating long after the fist packet has gone while 
the second has just arrived. 



3 Conclusions 

In this note we have considered neutrinos produced in the reaction e + A — > C + v. For 
the case of a large size of the electron wave packet (a e a a) we have calculated the 
neutrino wave packet (see eq.([16|)). Its size coincides with that of the incident electron 
wave packet. 

We have demonstrated that in the case of a e < p e \e — n\a a/Mq the suppression of 
neutrino oscillations at large distances is due to the neutrino energy dispersion, while in 
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the case of p e \e — n\aA/M c < cr e <C a a h occurs because of the neutrino packet separation. 
The corresponding suppression lengths are given by eqs.(J2Sl) and (J2*7j) . It is evident that 
for terrestrial conditions such lengths are unrealistic. 
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